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1 
Chapter   I 

INTRODUCTION 


The  rapid  growth  of  the  computer  industry  has 
provided  the  engineering  world  with  a  whole  array  of 
powerful  computing  machines  that  have  completely  altered 
most  aspects  of  engineering  practice  .  Computer  Aided 
Design  (CAD),  Computer  Aided  Manufacture  (CAM),  Computer 
Integrated  Engineering  (CIE)  and  Computer  Integrated 
Manufacture  (CIM)  are  some  of  the  technologies  that  are  a 
direct  result  of  the  computer  revolution. 

CAD  and  CAM  are  probably  the  fastest  expanding  areas 
in  industry  today.  These  two  fields  initially  emerged  as 
related  but  independent  technologies  .  A  large  volume  of 
current  research  is  focussed  on  efforts  to  link  these  two 
technologies  together  for  maximum  production  efficiency. 
Traditionally,  designers  have  worked  in  isolation,  taking 
only  design  considerations  and  performance  requirements 
into  account  while  creating  a  design;  manufacturing 
limitations  are  generally  not  given  much  consideration. 


Computer  Integrated  Manufacturing  (CIM)  is  a 
technology  developed  to  incorporate  CAD  and  CAM  into  a 
united  framework .  CIM  can  be  divided  into  two  categories 
based  on  the  function  it  performs  .  The  first  category 
deals  with  the  creation  and  manipulation  of  a  common 
database  that  can  be  used  by  the  designer  as  well  as  the 
production  planner.  This  is  a  very  significant  aspect  of 
CIM  and  plays  a  leading  role  in  streamlining  communica- 
tions between  the  design  and  manufacturing  departments 
within  a  company  .  Most  of  the  current  research  work  in 
the  field  of  CIM  falls  into  this  category  .  The  second 
category  is  generally  known  as  designing  for  producib- 
ility  .  The  motivation  here  is  the  existence  of  a  basic 
weakness  in  conventional  design:  even  though  the  design 
process  may  create  an  efficient  design,  it  does  not 
always  ensure  that  the  design  is  one  that  can  be  easily 
produced,  since  the  designer  is  operating  essentially 
independently  of  manufacturing  considerations.  In  conven- 
tional design  it  is  quite  common  to  sacrifice  performance 
to  obtain  better  producibility  .  However,  it  is  often 
possible  to  enhance  producibility  even  without 
sacrificing  performance  if  we  tackle  the  design  and 
manufacturing  problems  simultaneously.  This  approach. 


therefore,  exploits  the  idea  of  introducing  manufacturing 
considerations  at  the  design  level  itself. 

The  research  described  in  this  thesis  is  an  effort 
to  implement  the  "design  for  producibility"  aspect  of 
Computer  Integrated  Manufacturing.  The  objective  set  for 
this  work  was  to  devise  an  analytical  framework  for  the 
integration  of  CAD  and  CAM  and  to  develop  software  which 
allows  the  designer  to  implement  not  only  performance 
considerations,  but  also  manuf acturability  at  the  stage 
of  designing.  In  this  environment,  the  designer  includes 
manufacturing  considerations  such  as  manufacturing  time, 
manufacturing  cost,  production  sequences,  cutting  speed, 
cutter  diameter,  feed,  etc  .,  while  evaluating  and 
modifying  candidate  designs.  This  effort  will  broaden 
the  scope  of  the  design  problem  by  establishing  tighter 
links  between  the  design  and  manufacture  of  a  product;  it 
will  also  enable  designers  to  bridge  the  long-standing 
gap  between  performance  and  manufacturing  considerations. 

In  this  research,  the  specific  problem  of  designing 
and  manufacturing  a  plate  cam  is  considered  .  It  is 
assumed  that  the  cam  will  be  produced  by  end  milling  on  a 
Numerical  Control  Machine  .  The  specifications  of  a  Pratt 
&  Whitney  Tapemate  Series  C  NC  machine  were  taken  into 


consideration.  It  is  further  assumed  that  the  cam  will  be 
cut  from  a  rectangular  blank,  where  the  size  and 
orientation  of  the  blank  will  be  provided  by  the  program. 
Efforts  have  been  made  to  reduce  the  cost  of 
manufacturing  the  plate  cam  while  keeping  stringent 
constraints  on  the  performance  requirements. 

A  constrained  nonlinear  optimization  approach  has 
been  adopted  in  this  research  to  achieve  design/manufac- 
turing integration  .  This  approach  provides  a  convenient 
unified  framework  within  which  design  and  manufacturing 
considerations  can  be  dealt  with  simultaneously.  It  also 
provides  a  systematic,  iterative  method  for  choosing  the 
values  of  the  design/manufacturing  parameters  in  order  to 
create  a  product  that  is  optimal  with  respect  to  design 
as  well  as  manufacture. 

To  formulate  the  optimization  problem,  the  design 
engineer  is  required  to  specify  the  following: 

1)  A  set  of  design  variables  whose  values  must  be 
determined  to  obtain  the  optimum  solution  .  The 
engineer  is  also  required  to  specify  initial 


values  for  these  variables. 

2)  An   objective   function   to   be   minimized   or 
maximized. 

3)  Constraint  functions  that  represent  feasibility 
conditions  that  the  final  design  must  satisfy. 

The  type  of  cams  considered  in  this  research  are 
polynomial  cams,  i.e.,  cams  whose  profiles  are  specified 
by  a  polynomial  of  suitable  degree  .  In  conventional 
design,  the  number  of  coefficients  of  the  polynomials  are 
set  equal  to  the  number  of  boundary  conditions  specified. 
The  evaluation  of  these  coefficients  is  then  just  a 
matter  of  solving  a  set  of  linear  equations,  as  explained 
in  Chapter  2  .  The  values  of  the  coefficients  are 
therefore  uniquely  determined[ 1 ] ;  however,  the  design 
obtained  may  not  be  a  feasible  design  from  either  the 
manufacturing  or  the  performance  point  of  view  .  The 
approach  taken  in  this  research  is  to  choose  the  order  of 
the  polynomial  to  be  one  higher  than  that  required  for 
satisfying  the  specified  boundary  conditions.  This  extra 
coefficient  can  then  be  treated  as  a  design  variable  in 
the  optimization  problem. 


The  base  diameter  of  the  cam  plays  an  important  role 
in  cam-follower  mechanisms  .  The  minimum  possible  base 
diameter  is  usually  determined  by  the  size  of  the  roller 
follower,  the  pressure  angle  requirement  and  the  required 
lift  .  On  the  other  hand,  large  cams  are  undesirable 
because  of  the  bulkiness  of  the  mechanism  and  the 
presence  of  heavy  rotating  parts  .  The  base  diameter  is 
also  treated  as  a  design  variable;  however  the  upper  and 
lower  bounds  on  it  need  to  be  specified  .  As  an  option, 
the  base  diameter  can  be  treated  as  constant  during  the 
optimization  procedure,  if  the  user  wishes  to  do  so. 

After  establishing  the  design  variables,  the  second 
specification  in  a  nonlinear  optimization  technique  is 
to  define  a  cost  function  that  is  to  be  minimized.  Here, 
the  manufacturing  cost  is  considered  to  be  directly 
related  to  manufacturing  time;  hence  the  objective 
function  chosen  here  is  the  cutting  time  .  This  cost 
function,  however,  can  be  altered  as  desired  without 
significantly  changing  the  general  method  that  is 
developed  in  this  thesis. 

The  third  issue  in  the  optimization  process  is  the 
set  of  constraints.   Here,  both  the  designer  and 


manufacturer  can  impose  their  respective  constraints  on 
the  design  .  The  designer  may  constrain  the  design  by- 
putting  constraints  on  undercutting,  pressure  angle, 
velocity,  acceleration,  jerk,  contact  stresses,  etc, 
while  the  manufacturer  can  impose  constraints  on  cutting 
speed,  feed,  etc. 

Once  an  initial  set  of  values  is  given  for  the 
design  parameters,  the  coefficients  of  the  polynomials 
are  computed  using  the  boundary  conditions  and  the 
initial  values  assigned  to  the  design  variables  .  Now, 
using  these  polynomials,  the  constraints  are  evaluated  at 
specified  points  around  the  cam  contour.  The  tool  path  is 
then  generated  automatically  using  an  algorithm  that 
ensures  that  the  path  generated  will  be  the  most 
efficient  one  that  satisfies  the  manufacturing  tolerance 
requirements  .  After  the  path  has  been  generated,  the 
cutting  time  associated  with  the  path  can  be  estimated. 
As  mentioned  earlier,  this  cutting  time  is  the  cost 
function  .  If  all  constraints  are  satisfied  then  the 
estimated  cost  is  taken  directly  as  the  cost  function;  in 
case  there  are  violated  constraints  then  the  highest 
constraint  violation  is  multiplied  by  a  penalty  factor 
(usually  of  the  order  of  10  or  higher),  and  added  to  the 
estimated  cost  to  obtain  the  cost  function.   In  essence. 


if  any  constraint  is  found  to  be  violated  in  any  region, 
the  design  is  forced  to  move  away  from  that  region. 

Each  design  is  passed  through  a  grid  search  non- 
linear optimization  algorithm  in  order  to  compute  a  new, 
updated  design  .  By  cycling  through  a  suitable  number  of 
iterations,  the  design  variables  are  made  to  converge  to 
a  final  design  which  minimizes  the  cutting  time  and 
satisfies  the  constraints  that  are  imposed  .  The  final 
design  is  than  utilized  to  generate  Numerical  Control 
machine  codes  which  can  be  used  to  produce  the  cam. 

The  design  aspects  of  polynomial  plate  cams  are 
presented  in  Chapter  2  .  Manufacturing  aspects  and  the 
generation  of  the  tool  path  are  discussed  in  Chapter  3. 
Chapter  4  deals  with  the  formulation  of  the  nonlinear 
optimization  problem  with  emphasis  on  the  choice  of 
design  variables  and  the  formulation  of  cost/constraint 
functions  .  Some  of  the  numerical  examples  that  were  used 
to  test  the  proposed  model  are  presented  in  Chapter  5, 
and  some  conclusions  and  recommendations  for  future 
research  are  discussed  in  Chapter  6. 


Chapter  II 

PLATE  CAM  DESIGN 

Cams  are  used  to  lead  a  follower  element  by  direct 
contact  through  desired  motion  patterns  for  the 
performance  of  a  specific  task.  Cams  and  followers  may  be 
classified  in  different  ways  .  Based  on  physical 
structure,  cams  can  be  classified  as  wedge  cams,  barrel 
cams  or  plate  cams[l]  .  Similarly,  followers  can  be 
classified  as  roller  followers,  knife-edge  followers  or 
flat  face  followers.  This  work  focuses  on  plate  cams  with 
roller  followers  having  pure  translational  motions  . 
However,  with  proper  adjustments,  other  types  of  cams 
and/or  followers  may  also  be  incorporated. 

CAM  AND  FOLLOWER  MECHANISM 

Cam  and  follower  mechanisms  are  simple  and 
inexpensive;  they  have  few  moving  parts  and,  compared  to 
linkages,  they  occupy  less  space.  Cams  are  used  to  make 
the  follower  undergo  various  specified  motion  patterns, 
defined  in  the  form  of  rise,  return  and  dwell  phases, 
with  control  on  position,  velocity,  acceleration,  jerk. 
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etc  .  A  reciprocating  roller  follower  is  illustrated  in 
Figure  2.  The  follower  maintains  direct  contact  with  the 
cam  and,  as  the  cam  rotates  (generally  with  constant 
angular  velocity),  the  follower  moves  to  and  fro  along 
its  axis • 

The  displacement  diagram  (Figure  3)  provides  a 
clearer  understanding  of  the  performance  of  the  cam- 
follower  mechanism  .  Here,  the  abscissa  represents  one 
cycle  of  input  motion  (x)  .  The  ordinate  represents  the 
motion  of  the  follower  (y)  .  If  a  cam  is  rotating  with 
constant  angular  velocity  "w" ,  when  the  cam  rotates  by 
"x"  the  follower  will  move  a  distance  "y"  that  is 
dependent  on  the  cam  rotation,  i.e., 

y=f(x)  (1) 

The  overall  follower  motion  is  divided  into  four 
zones,  as  shown  in  Figure  3  .  The  region  where  the 
follower  is  moving  away  from  the  cam  center  of  rotation 
is  called  the  rise  zone.  The  maximum  rise  is  the  total 
lift  "L"  .  The  region  where  the  follower  is  approaching 
towards  the  cam  center  of  rotation  is  termed  the  return 
zone.  All  other  regions  are  dwell  zones.  The  dwell 


IS 


following  the  rise  is  called  the  top  dwell  and  the  one 
following  the  return  is  called  the  bottom  dwell  .  In  the 
dwell  zones,  the  follower  precisely  maintains  constant 
position. 

The  end  points  of  these  zones  are  defined  by  angles 
"  3-",  j=0,l,...,4,  as  shown  in  the  figure.  It  should, be 
ted  that    3,  is  the   starting  point  of  the  cycle,  and 


no 


0.  is  the  final  point  of  the  cycle;  hence,  the 
difference  between  them  should  always  be  equal  to  360 
degrees . 


With  these    6  . 's  known  "  aB-  ",  the  specific  range 


of  any  zone  i  can  be  computed  by; 


"i   '  'i    -    Bj-l  ''• 


The  procedure  by  which  a  relationship  is  developed 
between  the  follower  motion  and  cam  rotation  is  given 
below: 

The  displacement  of  the  follower  can  be  described  by 
an  equation  of  the  form. 

y=f(z)  (3) 
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Here,  "z"  is  the  dimensionless  ratio  of  the  angle  of 
rotation  of  the  cam  from  the  start  of  the  current  zone  to 
the  total  range  of  the  zone;  mathematically, 


z   =  e  /  Ag. 


(4) 


where 

9   is  the  angle  of  cam  rotation  measured 
from  the  start  of  zone  i. 

The  ratio  "z"  is  generally  termed  the  "time  ratio"  and 
varies  from  0  to  1 . 

It  is  assumed  that  the  follower  position  is 
normalized  by  dividing  the  true  position  by  the  total 
lift.  Thus,  the  normalized  position  "y"  in  Equation  (3) 
varies  from  0  to  1  .  Actual  position  can  therefore  be 
computed  by  multiplying  the  normalized  position  by  the 
total  lift  .  We  can  also  define  the  reduced  velocity  as 
the  derivative  of  position  with  respect  to  cam  rotation 
within  a  zone  "  9  ",  i.e., 

y  '  =  L  *  dy/d  9=L*f'(z)  (5) 
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where 

L  is  the  total  lift  of  the  follower. 

Similarly,  the  second  derivative  of  position  with 
respect  to  cam  rotation  within  a  zone  is  called  the 
reduced  acceleration  and  is  given  by  Equation  (6);  the 
third  derivative  of  position  with  respect  to  cam  rotation 
within  a  zone  is  called  the  reduced  jerk  and  is  given  by 
Equation  ( 7) . 


y"  =  L  *  d^y/d^e   =  L*f"(z)         (6) 


y"'=  L  *  d^y/d^0   =   L*f'"(z)         (7) 


The  true  position,  velocity,  acceleration  and  jerk 
of  the  follower  are  given  by  Equations  (8),  (9),  (10)  and 
(11)  respectively: 


y*   =  y  *  L  (8) 


where 


* 

y   is  the  true  position 

L   is  the  total  lift  of  the  follower 


y   =   y'   *  w  (9) 
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y    =    y"    *    w^  (10) 


y    =    y'"    *    w^  (11) 


The  radial  distance  from  the  cam  center  to  a  point 

on  the  cam  profile  is  given  by: 

r(z)=L*y  +  rj^  (12) 

where 

r(z)  is  the  radial  distance  to  the  cam  profile 

r,    is  the  base  circle  radius  of  the  cam 
b 

The  actual  follower  motion  depends  on  the  particular 
application  that  the  mechanism  is  designed  for  and  on  the 
boundary  conditions  specified.  It  is  important,  in  making 
the  transition  from  one  zone  to  next  ,  that  we  maintain 
continuity  of  position  and  its  derivatives  up  to  the 
degree  required  for  the  application  in  view.  In  order  to 
ensure  this  continuity,  boundary  conditions  are  enforced 
at  the  transition  points  between  adjoining  zones  .  The 
boundary  conditions  are  generally  expressed  in  terms  of 
normalized  position,  reduced  velocity,  reduced  accelera- 
tion, etc.,  and  are  specified  at  each  terminal  point  of 
each  rise  and  each  return  zone  .  For  example,  for  the 
follower  motion  shown  in  Figure  4,  the  boundary 


Figure   4:       Illustration  of  boundary   conditi 


ons 


20- 


conditions  are  specified  as  follows; 
RISE  ZONE 

z  =  0 

y  =  0,     y'=  0,     y"=  0 
z  =  1 

y  =  l,    y'=0,    y'"=0 


RETURN  ZONE 
z  =  0 

y=l,    y'=0,    y"'=0 
z  =  1 

y  =  0,     y'=0     y"=0 

There  are  various  techniques,  both  graphical  and 
analytical,  to  establish  a  relationship  between  the  input 
motion  of  the  cam  and  the  output  motion  of  the  follower. 
In  conventional  design,  standard  curves  are  adopted  to 
generate  the  required  cam  contour  for  specified  follower 
motions  and  boundary  conditions [ 2 ] .  Some  of  the  standard 
curves  that  are  commonly  used  are  cycloidal  curves, 
simple  harmonic  curves,  modified  harmonic  curves, 
parabolic  curves,  etc  .  The  criteria  for  selecting  the 
appropriate  curves  to  be  used  are  based  almost  entirely 
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Figure  5(a): 


Displacement  diagram  and  derivatives 
for  full-rise  modified  harmonic 
motion 

(Extracted  from  Reference  {1}) 


1    dl& 


Figure   5(b) 


Displacement  diagram  and  derivatives 
for  full-return  modified  harmonic 
motion  (Extracted  from  Reference   {1}) 
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Figure  5(c) 


Displacement  diagram  and  derivatives 
for  full-rise  cycloidal  motion 
:C£ttt«cted  from  Reference  {1}) 


Figure  5(d) 


Displacement  diagram  and  derivatives 
for  full-return  cycloidal  motion 
(Extracted  from  Reference  {1}) 
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Figure   5(e): 


Displacement  diagram  and  derivatives 
for  full-rise  simple  harmonic 
motion 

(Extracted  from  Reference  {1}) 


Figure   5(f) 


Displacement  diagram  and  derivatives 
for  full-return  simple  harmonic 
motion 

(_Extracted  from  Reference  {1}} 
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on  the  boundary  conditions  to  be  satisfied.  Figures  5a-5f 
illustrate  the  displacement  diagrams  and  derivatives  for 
full-rise  and  full-return  for  some  of  these  standard 
motions[l]  .  If  the  follower  is  required  to  have  zero 
acceleration  at  each  terminal  point  of  each  zone  then 
cycloidal  motion  will  be  the  obvious  choice,  whereas  if 
the  jerk  is  required  to  be  zero  at  each  terminal, 
modified  harmonic  motion  will  be  selected  .  This  is  a 
tedious  procedure  since  each  different  case  must  be 
implemented  by  a  different  curve.  Moreover,  apart  from 
being  cumbersome,  this  procedure  restricts  the  range  of 
possible  solutions  since  standard  curves  are  only  a  small 
subset  of  the  set  of  possible  solutions  which  satisfy  the 
specified  boundary  conditions. 

POLYNOMIAL  CAMS 

The  approach  adopted  in  this  research  is  to  use 
polynomials  to  define  the  cam  profile  .  Polynomials 
provide  a  convenient  means  by  which  a  wide  variety  of 
motion  patterns  can  be  generated[3]  .  Even  complicated 
cam  profiles  which  would  be  difficult  to  generate  using 
standard  curves  can  be  conveniently  described  using 
polynomials  .  By  merely  adding  successive  higher  order 
terms  to  the  polynomial,  different  curves  can  be 
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generated  .  These  polynomials  can  be  employed  as  single 
continuous  functions,  or  may  comprise  parts  of  a 
piecewise  defined  function[3].  In  the  current  work,  it  is 
assumed  that  two  different  polynomials  are  used  to 
represent  the  rise  and  return  respectively. 

The  polynomials  considered  in  this  work  are 
functions  of  the  time  ratio  "z"  .  Any  number  of  boundary 
conditions,  specified  with  respect  to  z,  can  be  satisfied 
by  using  an  appropriate  order  of  polynomial  .  The  order 
of  the  polynomial  is  generally  taken  to  be  n-1,  where  n 
is  the  total  number  of  boundary  conditions  specified. 

A  typical  rise  or  return  polynomial  can  be  written  in  the 
form: 


y(z)=c^+C2*z+C2*z^+c^*z-^+ +c  *z^"~-'-^     (13) 


The  derivative  of  this  polynomial  with  respect  to- 
"  0  "will  give  the  reduced  velocity: 


y'{z)=L/Ag  *(C2+2*C2*z+3*c^*z^+~+(n-l)*c  *z^""^^ )   (14) 
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where 

L   is  the  total  lift  of  the  follower. 

Similarly,  the  second  derivative  of  the  polynomial 
will  yield  the  reduced  acceleration,  and  the  third 
derivative  will  give  us  the  reduced  jerk: 


y"(z)=L/  Ae*{2*C3+6*C4*z+— +(n-l)*(n-2)*c^*z^"  ^^  )       (15) 
y'"(z)=L/  Aa*(6*c^+— +(n-l)*(n-2)*(n-3)*c^*z^""^^)  (16) 


Equations  (13),  (14),  (15)  and  (16)  can  be  used  to 
compute  the  coefficients  of  the  polynomial  once  the 
boundary  conditions  are  known. 

The  radial  distance  to  the  cam  profile  "r"  can  now 
be  expressed  by: 


r(z)=L*(cT+c„*z+c,*z^+c  *z^+ +c  *z''^~^^  }+r,  (17) 

i   z    J     4  n  D 


where 


r,   is  the  base  circle  radius  of  the  cam 


Note  that  the  first,  second,  and  third  derivative  this 
radial  distance  are  equal  to  y',  y'',  y'''  respectively. 
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EVALUATION  OF  POLYNOMIAL  COEFFICIENTS 

In  conventional  cam  design  the  number  of  free 
coefficients  in  the  polynomial  is  taken  to  be  equal  to 
the  number  of  boundary  conditions  specified.  In  this  case 
the  coefficients  are  uniquely  determined  through  the 
solution  of  a  set  of  simultaneous  linear  equations[ 1 ] .  A 
modified  version  of  this  approach  is  utilized  in  this 
work  and  is  explained  below. 

Instead  of  setting  the  order  of  the  polynomial  to  be 
equal  to  be  one  less  than  the  number  of  boundary 
conditions  it  must  satisfy,  we  choose  it  to  be  of  the 
same  order  as  the  number  of  boundary  conditions  .  The 
leading  coefficient  can  then  be  chosen  arbitrarily  and 
can  therefore  be  treated  as  a  design  variable  .  The 
solution  scheme  is  setup  in  such  a  way  that  any  change  in 
this  design  variable  will  automatically  effect  a  change 
in  the  remaining  coefficients  of  the  polynomial  so  that 
the  boundary  conditions  are  enforced  .  The  process  by 
which  this  is  achieved  is  the  following: 

1)  Define  the  boundary  conditions. 

For  example,  let  us  consider  the  boundary  conditions 

V 
for  a  rise  cycle  only: 
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At  z  =  0: 

y  =  0,     y'=  0,     y"=0. 
At  z  =  1: 

y  =  l,    y'=0,    y"'=0. 
2)  Select  a  polynomial  of  order  one  less  than  the  number 
of  boundary  conditions: 


2      3      4      5 
y=c, +c„*z+c-,*z  +c.*z  +Cc.*z  +c,*z 
-'123      4      5      5 

3)  Raise  the  order  of  this  polynomial  by  one: 

y=c, +c_*z+c,*z  +c.*z  +c^*z  +Cc*z  +c^*z 

4 )  Compute  y',  y''  &  y''': 

y'=C2+2*C3*z+3*c^*z^+4*C5*z^+5*Cg*z'^+6*c^*z^ 

y"=2*c^+6*c  *z+12*Cc*z^+20*c,*z^+30*Ct*z^ 
3     4         5  6  7 

y"  '=6*c^+24*C5*z-i-60*Cg*z^+120*c^*z^ 


(18) 

(19) 

(20) 
(21) 
(22) 


5)  Write  the  boundary  conditions  in  matrix  form,  using 
the  above  Equations  to  substitute  for  y,  y',  y''  and 

^    ^    * 

Y    : 


1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

2 

3 

4 

5 

6 

0 

0 

0 

6 

24 

60 

120 

-^7- 


0 
0 
0 

1 

0 
0 


(23) 
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In  general,  this  matrix  equation  can  be  simply  written 
as : 


I  I    ^1      I  I    qi    I 

A     I  I   :     I    =  I   :    I         (24) 

-I  I-  S.l-I  I-  ^n  J 
n,n+l           n+l  n 


The  matrix  A  can  be  partitioned  as 


(25) 


n,n+l  n,n+l 

where 

[V]  is  an  (n  X  1)  column  vector  that  is 

identical  to  the  (n+l)th  column  of  [A] 
[TT]  is  the  (n  x  n)  submatrix  consisting  of  the 
first  n  columns  of  [A] 
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6)  Now,  Equation  (24)  can  be  manipulated  into  the 
following  form: 


n  — ' 


n,n 


I  +  c^+i  I   V   I   (26) 


'n  — ' 


n 


n 


n,l 


7)  The  last  step  is  to  solve  this  system  of  linear 

equations  to  obtain  c,  through  c, . 

i  6 


1 

0 

0 

0 

0 

0 

^1 

0 

0 

0 

1 

0 

0 

0 

0 

^2 

0 

0 

0 

0 

2 

0 

0 

0 

C3     1  = 

0 

^    ^7 

0 

1 

1 

1 

1 

1 

1 

^4 

1 

-1 

0 

1 

2 

3 

4 

5 

1    C3 

0 

-6         1 

_0 

0 

0 

6 

24 

60_ 

-^6- 

0_ 

_-120_ 

For  example,  let  the  design  variable  be  equal  to  1 


i.e. ,   c^  =  1 


Then,  after  solving  the  above  system  of  equations,  the 
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final  polynomial  coefficients  obtained  are: 


^1  =  ^2  =  ^3  =  °- 

c^  =  5.3326 

C(.  =  -4  .  666 

c,  =  -.6666 
6 


DESIGN  CONSIDERATIONS 

Before  setting  up  the  design  process,  we  must  first 
establish  a  set  of  important  design  considerations. 
In  the  current  work,  these  design  considerations  are 
enforced  through  constraints  in  the  optimization  process. 
Design  considerations  are  user  dependent  and  may  vary 
from  case  to  case  .  Some  basic  design  considerations  are 
explained  below;  however,  the  user  may  add  performance 
criteria  as  desired. 

Pressure  Angle 

One  of  the  most  critical  parameters  in  cam  design  is 
the  pressure  angle  .  The  pressure  angle  is  defined  as  the 
angle  between  the  axis  of  the  follower  stem  and  the  line 
of  action  of  the  force  exerted  by  the  cam  on  the  roller 
follower.  High  values  of  pressure  angle  should  be  avoided 
in  order  to  reduce  friction,  minimize  chattering  of  the 
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Figure   6 : 


Plate  Cam  with  reciprocating  roller 

f ol lower 

(Extracted  from  Reference  {17}) 
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follower,  and  to  obtain  better  force  transmission[ 2 ] .  As 
can  be  seen  in  Figure  6,  only  the  force  component  along 
the  axis  of  the  follower  stem  is  used  to  overcome  output 
load.  By  keeping  the  value  of  $  as  low  as  possible  the 
tangential  component  (unused  component)  of  force  is 
reduced,  thus  reducing  the  sliding  friction  between  the 
follower  and  its  guideway  .  The  distance  from  the  cam 
center  of  rotation  to  the  trace  point,  "Ho"  is  given 
by[4]: 

Hq  =  sqrt((rj^  +  r^)^-e'^)  (27) 

where 

r,  is  the  base  circle  radius  of  the  cam 
b 

Tr:    is  the  roller  follower  radius 
e   is  the  follower  offset 

The  instantaneous  distance  between  the  trace  point  and 
the  cam  center  of  rotation,  "H"  is  given  by, [4] 

H  =  Hq  +  y  *  L  (28) 

where 

y  is  the  normalized  position  of  the  follower 

L  is  the  total  lift  of  the  follower 
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Now,  the  pressure  angle  is  given  by: 

*   =  tan"-^[(y'-e)/H]  (29) 

where 

y'  is  the  reduced  follower  velocity  at  any 
instant 

The  maximum  pressure  angle  should  generally  be  less 
than  30  degrees  for  most  applications. 

Undercutting 

The  curvature  of  the  cam  is  another  important 
consideration  in  cam  design .  If  the  curvature  of  the  cam 
is  too  high,  the  follower  may  not  be  able  to  track  the 
surface  precisely [ 1 ] ;  there  will  also  be  an  increase  in 
the  contact  stresses  developed.  This  phenomenon  is  called 
undercutting  and  must  be  avoided  if  the  cam  is  to  perform 
satisfactorily.  With  a  roller  follower,  the  undercutting 
occurs  in  a  convex  curve  when  the  radius  of  the  roller 
follower  is  greater  than  the  radius  of  curvature  of  the 
cam  profile[2] .  With  a  concave  curve,  it  occurs  either 
when  the  radius  of  curvature  of  the  pitch  curve  is  less 
than  zero  or  when  it  is  less  than  the  radius  of  the 
roller  follower.   For  a  safe  design,  the  radius  of 
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Figure   7 ; 


Illustration  of  Undercutting 
(Extracted  from  Reference  {16}) 
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curvature  of  the  profile  must  be  greater  than  one  and  a 
half  times  the  radius  of  the  roller  follower[4]. 

The  radius  of  curvature  of  the  cam  is  given  by: 


[H^  +  (y'-  e)2]^/2 


H^  +  (y'  -  e) (2*y'  -  e)  -  H 

o 


(30 


There  are  two  different  ways  to  avoid  undercutting[ 2 ] : 

1)  Use  a  roller  follower  with  small  radius.  In  this 
case,  however  there  will  be  high  stresses  in  the 
follower. 

2)  Use  a  large  base  radius  for  the  cam.  However,  if 
the  cam  is  too  large  it  will  create  vibration 
problems  and  increase  space  requirements. 

In  this  work  the  roller  follower  radius  is  kept 
constant  and  therefore  undercutting  effectively  imposes 
limits  on  the  base  radius. 
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Separation 

The  above  mentioned  pressure  angle  and  undercutting 
considerations  are  very  important  in  the  design  process 
as  they  impose  limits  on  cam  size  .  Another  essential 
consideration  in  the  design  of  a  cam-follower  mechanism 
is  the  assurance  of  continuous  contact  between  cam  and 
follower  surfaces  at  all  times  during  the  cycle  -  Here, 
we  are  required  to  take  a  look  at  the  dynamics  of  the 
system  (Figure  6).  The  driving  force  "F,"  is  given 
by: 


F^  =  (  Pj^  +  m^*Y    )/cos  $        (31) 
where 

P,  is  the  preload  specified  by  the  user 

(including  friction,  spring  force,etc.). 
m^  is  the  mass  of  the  follower  mechanism, 
y   is  the  true  acceleration  of  the  follower. 


Physically,  the  force  F,  can  only  act  upwards  on  the 
follower.  If  the  follower  acceleration  required  to  track 
the  cam  profile  on  the  return  is  too  high,  a  negative 
value  of  F,  might  be  required  to  maintain  contact.  Since 
this  is  not  possible,  this  situation  will  result  in 
separation  between  the  follower  and  cam  surfaces,  and 
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this  is  generally  unacceptable  .  In  order  to  avoid 
separation  the  resultant  driving  force  is  required  to  be 
positive  all  through  the  cycle  .  One  way  to  prevent  the 
driving  force  from  going  negative  is  to  increase  the 
preload[6],  but  that  will  result  in  an  increase  in  the 
contact  stress  developed  .  The  best  way  to  avoid 
separation  is  to  keep  the  acceleration  on  the  return  as 
low  as  possible  by  controlling  the  profile.  In  this  work, 
separation  between  the  follower  and  cam  surfaces  is 
avoided  by  restricting  the  driving  force  to  be  positive 
throughout  the  motion . 

Stress  Considerations 

If  we  consider  an  ideal  case  of  two  convex  rigid 
bodies,  we  will  have  a  case  of  straight  line  contact. 
However,  in  practice  the  bodies  deform  to  give  area 
contact.  As  studied  by  Hertz,  stress  distribution  between 
two  contacting  bodies  with  cylindrical  shape  and  with 
perfect  alignment  is  such  that  the  compressive  stress  is 
a  maximum  at  the  center  of  the  contact  area  and  decreases 
to  zero  at  the  ends;  however,  in  actuality,  deflection 
and  misalignment  may  lead  to  extremely  high  local  stress 
generation[2 ] .  In  this  work,  the  contact  stress  is 
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computed  by  using  the  Hertzian  Equation[4] 


F   (  y  +   r   ) 
2 
^C     =  (32) 


Stress  is  a  very  critical  factor  because  many 
problems  of  high  wear  and  failure  are  caused  by  excessive 
stresses;  consequently,  the  designer  is  required  to  take 
steps  to  minimize  the  stresses  encountered. 

Elimination  of  oscillation  in  the  follower  velocity 

The  use  of  high  order  polynomials  tends  to  produce 
oscillations  in  the  path.  These  oscillations  are  highly 
undesirable  in  cam  follower  mechanisms.  In  order  to  avoid 
these  fluctuations,  constraints  are  imposed  on  the 
velocity  of  the  follower  to  ensure  a  steady  rise  and  a 
steady  return  .  During  the  rise  cycle  it  is  required  that 
the  follower  velocity  always  be  positive,  whereas 
during  the  return  cycle  it  should  always  be  negative. 

All  the  above-mentioned  design  considerations  are 
enforced  as  design  constraints  which  need  to  be  satisfied 
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by  the  final  design  .  However,  the  user  may  include  any 
other  design  considerations  that  he  wishes  to  use  as 
constraints  in  the  optimization  problem  .  The  exact 
constraint  equations  that  are  used  to  implement  the 
considerations  discussed  above  are  outlined  in 
Chapter  4 . 
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Chapter  III 


MANUFACTURING  ASPECTS 


The  production  of  cams  using  NC  machines  was  started 
in  the  mid  50 's.  Since  then,  digital  computers  have  been 
utilized  extensively  to  design  cams  and  to  generate  data 
points  for  machining  the  cam  profile.  These  data  points 
are  transferred  onto  a  metallic  tape  or  other  data 
transfer  medium,  which  is  then  fed  to  the  NC  machine  . 
Here,  except  for  the  required  machine  setup,  human 
involvement  is  almost  completely  eliminated[ 7 ]  .  The 
automation  of  cam  production  has  not  only  reduced  the 
manufacturing  time  drastically,  but  has  also  increased 
the  quality  of  the  product. 

PRODUCTION  PROCESS 

For  the  purpose  of  this  research,  the  NC  machine 
considered  is  a  Pratt  &  Whitney  Tapemate  Series  C 
Numerical  Control  machine  .  This  machine  can  perform 
operations  like  milling,  drilling,  boring  and  tapping 
economically  and  ef f iciently [ 8 ]  .  Table  1  provides  the 
relevant  specifications  of  this  machine.  The  end  milling 
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Table      1 


Numerical    Control    Machine 
Specifications 


Mclrlr    f:yKlcn 


Inch    Sy?;I  cm 


Ciip.Tclty 


Dr II  line 


Tapping 
llthinR" 


rc;5    C.l5t    Iron 
S4SC  Mild   steel" 


_KC25_C.TS£  Iron  _ 
SisCMlld   st:cr:i_ 
"tjid   mtil 


M.'i>i.     20iiini  ili»._ 

3  "^    Ifiinm  (J  Ui.  _ 
3  '^    1  Otnin   tl  In. 
M.1X.    iSnim  dla. 


tl.ix.     i/'.''  .li.i. 
Il.i>:.     l/'i"   di.i. 
1/S  ^    Vii""   .li.iV 
1/tl  •>     '■'/li';  di:'- 
MctX.     5/;!"    dia. 


Tabic 


_ Servo   nolor 

Sl7i- 

T''L*J  '^_^T.'?*'^J- 

l.'j.Ttt    coimnnnd    Lncrewienc 

Rapttl    traverse 

Feed    I  ate 


I'osit ionlnp.  accuracy 
P-epeatal>llity 


Max.    weight  on    tabic 
Splndic  head   T^- axis}"" 
Servo  motor 


FANliCJlCsi-rv^O.  motor    mod 
fiin'x    '.f-O    mill 
')(")   >:_.  isd  ifm 

IJ.Ul   mill/ pill v.e_ 

6000   mm/min.      

r '-    2000   limi/itiLn.         


I   .'. 


.1/1 


in 


I ' 


?/«:•_ 


iO.05    mmndOjiim 
iO.'ol'nm 


150   kc 


ii.iiiiii,i'7pulsc 

2S()     i|.m 

ll.O:>  ■      100    ipm 

■n.'(iir'"/i2''    ^ 

io.liUUj'^ 

J  j()     Ui5 


Spindle    travel  

nlstancc    from   spindle   nose 
t-o   table   surface 


Least   comniand   Increment 

I^apid    traverse 

Feed^atr      ^     

^J^sii  tinning  accuracy 

l'.e|ieat.Tli^Htj  

tJunber   o(    tools 


FANUC  I'C   scrvii  mo^or  moilei    ij 
2  20  mm 


130  •^    350 


5  •'•   11.5" 


(^nnttnVi"d    ■*t>:i*n 
Slmii  I  f  aiiooii*:  I  y 
__r(in_t  ipl l;ib  1  o_  axes 

Intfri'Olat  lt.in_ 

1-ixcd    cvclc 


Tool    lpnj;rh   offset 

Aool   radius  offset _ 

J^Veir"'^ . '. 

_IUif  fcr_5t_oraj;q 

Manual    rt'Iurn    to    rcro 

__Auromatic    rrturn  _to   ?.cto_ 
_nt]}    anii_»n.        " 
_Total    woip.ht  (Approx.  ) 

Power  supply 


Basic  Option 


Spindle 
Motor 


FANUC  DC  spindle 
motor  model  2 


Pole  change  type 
induction  motor 


_0_^1   mm/^lso 

I'lOOO   nim/min .      

"l   ^\^ilOO  'iW'/riin._ 

'Pj!!"^'  iT^/JCiyjra _ 

•_0_.01_mm 

_7_J'_^ L    ."__-     '-" 

i,    \,  }.      '  '    .. 

2    axes    of    X   >    or    X-7. 

I. lunar    and    (1  f  rriilar_  i  nf  ■■  r  |- 
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slotting  process  is  selected  for  manufacturing  the  cam. 
It  is  assumed  that  the  cam  will  be  cut  from  a  rectangular 
blank,  the  dimensions  of  which  will  be  provided  by  the 
optimization  process. 

In  practice,  the  selection  of  milling  cutters  for  a 
certain  task  is  totally  based  on  the  engineer's 
experience.  For  this  work,  the  required  machining  data 
was  acquired  from  [9]  .  This  reference  provides  a  table 
which  prescribes  preferred  values  of  machining  parame- 
ters for  end  mill  slotting  on  alloy  steel  material  . 
According  to  this  table,  a  tool  diameter  of  3/4"  is  the 
most  appropriate  from  the  point  of  view  of  both  tool  life 
and  productivity.  This  is  also  evident  from  table  2,  as 
the  3/4"  diameter  tool  provides  a  good  feed  rate  and 
cutting  speed  for  a  wide  range  of  materials.  Some  of  the 
prescribed  machining  data  obtained  from  [9]  is  given 
be 1 ow : 

Tool  material High  Speed  Steel 

Tool  diameter     3/4   inches 

Number  of  teeth   4 

Depth  of  cut     0.25  inches 

Feed  per  tooth    0.002inches 

Cutting  speed    110  feet  per  minute 
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Once  the  cutting  speed  is  set  and  the  tool  diameter 
is  known,  the  required  spindle  speed  "V  "  in  rpm  can  be 
calculated  from[10] 

V  =(12*V  )/(^  *  d^)  (33) 

where 

V  is  the  nominal  cutting  speed  in  fpm 
d.  is  the  tool  diameter  in  inches 
There  are  only  six  fixed  spindle  speed  available  on 
the  NC  machine,  so  the  closest  lower  speed  (V  ')  is 
automatically  selected.   Using  this  spindle  speed  actual 
cutting  speed  can  be  computed  from 


V  =12/(60*  ""     *d^)  *  V^'  (34) 

ca  t     s 

where 

V  is  the  actual  cutting  speed  in  fpm 

V  '  is  the  actual  spindle  speed  in  rpm 


Also,  the  feed  in  inches  per  second  can  be  computed 
from  the  known  spindle  speed[14]: 

Feed  =  fpt  *  N^  *  V  V60  (35) 

t    s 

where 

fpt  is  feed  per  tooth  in  inches 

N.   is  the  number  of  teeth  in  the  tool 
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TURNS  IN  TOOL  PATH 

In  computing  the  cutting  time  it  was  anticipated 
that  the  number  of  turns  in  the  tool  path  would  be  a 
significant  factor.  A  simple  test  was  conducted  to  verify 
this  hypothesis  and  to  obtain  data  regarding  the  time 
lost  in  these  turns.  The  details  of  this  test  are  given 
be 1 ow : 

The  basic  idea  of  the  test  was  to  measure  the  time 
required  for  the  tool  to  travel  a  fixed  distance  through 
different  paths  .  All  other  parameters  such  as  the  feed 
rate  and  speed  are  kept  constant  .  The  NC  machine  was 
first  programmed  so  that  the  tool  traversed  the  sides  of 
a  square  of  side  4  inches.  The  machine  was  then  programm- 
ed to  drive  the  tool  from  corner  1  to  corner  3  and  back 
to  corner  1  along  the  staircase  path  shown  in  Figure  6. 
This  path  has  20  turns  in  the  tool  path  .  Similar 
staircase  paths  were  generated  with  40,  80  and  100  turns. 
Note  that  the  total  length  of  all  these  paths  is  equal. 

In  each  case,  the  tool  was  driven  along  the 
generated  path  for  20  cycles  .  The  results  obtained  are 
shown  in  Table  3.  Based  on  these  results,  an  average  time 
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Figure   8 : 


Estimation  of  time  lapse  per  turn  in 
tool  path 


Table   3 ; 


48 

Estimation  of  Time  Lapse  per  turn 


Case   1 


No.  of  cycles 


No.  of  turns/cycl    Total  time  (Min) 


20 
20 


4 
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14.155 
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5.216 
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25.58 


.3  3seconds 
.  39seconds 
.  37seconds 
.  38seconds 
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lapse  of   0  .4   seconds  per  turn  was  assumed   in  all 
subsequent  calculations. 

GNERATION  OF  DATA  POINTS 

The  generation  of  data  points  for  the  tool  path  for 
manufacturing  is  a  critical  and  time  consuming  step  in 
process  automation  .  Conventionally,  a  series  of  fixed 
steps  is  taken  along  both  the  x-axis  and  the  y-axis  to 
generate  the  cam  contour.  A  detailed  description  of  this 
method  is  presented  in  Reference  [6]  .  The  size  of  this 
fixed  step  is  determined  at  the  most  critical  section  of 
the  cam  profile,  where  the  step  size  is  usually  required 
to  be  very  small  in  order  to  meet  the  desired  accuracy. 
This  approach  has  been  modified  in  this  work  .  It  is 
assumed  that  it  is  not  necessary  to  maintain  the  same 
step  size  throughout  the  cam  contour  because  at  some 
sections  of  the  cam  a  larger  step  size  may  be  permissible 
for  the  same  accuracy  .  The  advantage  of  this  variable- 
step  strategy  is  that  it  minimizes  the  number  of  turns  in 
the  tool  path,  thus  reducing  the  manufacturing  time 
significantly. 

A  Newton-Raphson  iterative  procedure  is  adopted  to 
generate  the  variable  steps  along  the  x  and  y  axes  over 
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Figure      9 : 


Fixed  step  path  generation 
(Extracted  from  Reference  {6}) 


Figxire  10: 


Variable  step  path  generation 
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the  cam  contour  .  It  is  assumed  that  the  movement  of  the 
tool  is  restricted  to  be  on  the  outside  of  the  cam 
profile  at  all  points  .  In  other  words  the  cutter  never 
cuts  into  the  cam  profile,  but  remain  outside  the  actual 
profile  within  a  pre-specif ied  accuracy  "  e",  as  shown  in 
Figures  11a  and  lib.  This  strategy  makes  it  possible  to 
improve  the  accuracy  further  by  later  passing  the  cam 
through  finishing  processes,  such  as  profile  grinding. 
The  variable  step  generation  technique  developed  is 
described  below: 

1)  The  entire  profile  of  the  cam  is  divided  into 

zones  based  on  the  sign  of  the  slope,  as  shown  in  Figure 
10  .  This  is  done  by  computing  the  critical  points  with 
respect  to  the  coordinate  axis,  i.e.  points  at  which  the 
slope  is  either  zero  or  infinity  .  Between  any  two 
consecutive  critical  points,  the  slope  "dy/dx"  maintains 
the  same  sign. 

2)  When  all  the  critical  points  are  known,  the  process 
of  data  point  generation  begins  at  the  first  critical 
point  and  follows  the  cam  contour  until  it  comes  around 
to  the  same  starting  point.  As  explained  above,  the  tool 
path  is  required  to  lie  outside  of  the  cam  contour.   For 
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example,  if  slope  is  positive,  the  tool  should  first  take 
a  step  away  from  the  cam  contour  along  the  x-axis  and 
then  step  towards  the  cam  contour  along  the  y-axis.  The 
radial  distance  between  the  data  point  and  the  cam 
profile  must  be  equal  to  the  required  accuracy  "£". 

This  approach  is  implemented  as  follows: 

If  dy/dx  >  0 

f(0)  =  (r(0)  +  e)  *  sin(0)   -  y  =  o         (36) 
g{a)   =   r(a)   *  cos(a)   -  x,  =  0  (37) 

Else  if  dy/dx  <  0 

f(0)  =  (r(Q)  +  e)  *  cos(0!)j  -  x  =  0         (38) 
g(a)  =   r(a)   *  sin(a)  -  y^^  =  0  (39) 

where 

e  and  a   are  the  angle  of  cam  rotation 

rf0)  is  the  instantaneous  radial  distance  to  the 
cam  profile 

In  order  to  generate  the  tool  path,  we  only  need  to 
solve  the  above  equation  for  0  and  o  at  each  step  . 
Referring  to  Figure  11(a),  the  slope  in  section  1  is 
negative.  Suppose,  x^  and  y^  are  known  coordinates  of  a 
point  on  the  profile,  and  the  angle  is  known.  Starting 
from  this  position  we  want  to  move  along  the  y-axis  to 
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Figure  11(a):    Generation  of  data  points 
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Figure  11(b):    Generation  of  data  points 
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(x,  ,y,  )  such  that  the  maximum  radial  deviation  from  the 
true  profile  is  "  g  "  .  The  new  angle  9  is  computed  by 
solving  Equation  38  using  an  iterative  Newton-Raphson 
procedure  described  by 

^+1  =%  -  f(  ^)/f'<  %)         (40^ 
where 

n  is  the  iteration  number 

and  f '(  0  )  is  given  by; 


f'(  ^^)=r'(   e^)  *  cos(  0j^)-(r(  0j^)  +  e  )  *  sin(  e^)   (41) 


By  using  this  iterative  procedure  the  angle  "  e  "  can 
be  obtained  .  This  will  enable  us  to  calculate  the  step 
along  the  y-axis  .  Next  the  tool  must  be  moved  to  the 
position  (x2/y2)  shown  in  the  figure.  This  is  achieved 
by  again  using  a  Newton  Raphson  iteration  to  find  the 
angle  " a  "  in  Eq.  39: 

«n+l  =    °n  -  5^  "n^/^'^  "n^        ^^2) 


g  '("n^  "  ^'^  "n^  *  sin(  a  ^)  +  r (  o^ )  *  cos(a^)  (43) 


In  the  same  manner,  data  points  can  be  generated  in 
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sections  of  the  cam  contour  where  the  slope  is  positive 
by  solving  from  the  functions  given  in  Equations  36  and 
37. 

The  total  number  of  steps  generated  is  recorded 
along  with  the  data  points  .  The  total  number  of  turns  in 
a  tool  path  is  considered  to  be  equal  to  the  total  number 
of  steps  taken  by  the  tool. 

THE  ESTIMATION  OF  CUTTING  TIME 

Once  the  data  points  for  the  tool  path  have  been 

generated  for  the  entire  cam  contour,  the  distance  "S.  " 

■cp 

that  the  tool  will  travel  in  order  to  cut  the  cam,  can  be 
computed  as  follows: 

Stp=    t<-j  -  ^j-l)'  -^  (^j  -  ^j-l^']'^'     ^''^ 


where 


X.  is  the  x-coordinate  of  data  point  "j" 
y.  is  the  y-coordinate  of  data  point  "j" 


When  the  feed,  total  distance  travelled  and  total 
number  of  turns  are  known,  the  cutting  time  can  be 
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estimated  as  follows: 


F  =   S^  /Feed   +  0.4  *  N„.„  (45) 

tp'  stp 


where 


N  .   is  the  total  number  of  turns, 
stp 


This  estimate  of  cutting  time  is  then  used  as  the 
cost  function  in  the  optimization  problem  that  will  be 
discussed  in  detail  in  the  next  chapter. 
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Chapter  IV 

OPTIMIZATION 

Nonlinear  optimization  techniques  provide  a 
convenient  means  through  which  design  objectives  can  be 
achieved;  by  using  suitable  constraints,  we  can  be 
assured  that  the  design  satisfies  the  required  perform- 
ance criteria  .  In  order  to  establish  an  optimization 
problem,  the  optimality  criteria  need  to  be  defined.  In 
order  to  do  this,  the  following  are  required  to  be 
identified: 

1)  A  set  of  design  variables 

2)  A  cost  function  that  is  to  be  minimized 

3)  A  set  of  constraints  that  must  be  satisfied 

1)  CHOICE  OF  DESIGN  VARIABLES 

The  first  and  foremost  problem  is  to  pick  out  the 
right  set  of  design  variables  to  work  with  .  A  thorough 
investigation  is  needed  to  determine  which  parameters 
should  be  treated  as  design  variables  and  which 
parameters  should  be  kept  constant.  It  may  appear  safe  to 
treat  all  parameters  in  sight  as  design  variables,  but 
for  diverse  reasons  this  is  not  possible  or  practical [ 11 ] 
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Material  selection,  for  instance,  is  not  treated  as 
design  variable  in  practice;  rather  it  is  selected  on  the 
basis  of  experience  by  the  engineer.  The  number  of  design 
variables  must  also  be  reduced  as  much  as  possible  in 
order  to  avoid  computational  complexities  and  to  cut  down 
computing  time.  In  this  research,  quite  a  few  parameters 
could  have  been  easily  treated  as  design  variables,  for 
example,  tool  diameter,  cutting  speed,  roller  diameter, 
cam  and  roller  thickness,  etc  .  Doing  so  would  have 
resulted  in  a  much  more  complex  problem,  which  would  have 
been  very  expensive  to  optimize;  in  certain  cases, 
reaching  the  optimum  solution  may  not  even  have  been 
possible . 

The  approach  taken  in  this  work  was  to  treat  the  cam 
diameter  and  the  orientation  of  the  cam  on  the  NC  machine 
table  as  design  variables  in  addition  to  the  extra  coeff- 
icients in  the  rise  and  return  polynomials  .  This 
combination  of  design  variables  provides  a  convenient  and 
practicable  approach  to  optimization  in  the 
design/manufacture  environment  .  These  design  variables 
are  explained  in  detail  below. 
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The  leading  polynomial  coefficients: 

The  polynomial  is  a  convenient  tool  to  describe  a 
set  of  points  in  a  plane  .  The  order  of  the  polynomial 
depends  on  the  number  of  conditions  it  is  required  to 
satisfy,  i.e.,  the  number  of  boundary  conditions  in  our 
case  .  As  explained  in  Chapter  2,  the  approach  taken  in 
this  work  is  to  choose  the  order  of  the  rise  and  return 
polynomials  to  be  one  higher  than  the  order  required  to 
satisfy  all  the  boundary  conditions  .  The  leading  coeff- 
icients of  these  polynomials,  which  are  treated  as  design 
variables,  can  be  chosen  arbitrarily;  the  remaining 
coefficients  can  then  be  calculated  in  order  to  satisfy 
the  boundary  conditions.  The  optimization  process  will 
choose  values  of  the  leading  coefficients  in  such  a  way 
that  feasibility  and  cost  minimization  are  achieved  . 
These  design  variables  usually  have  a  strong  effect  when 
the  design  is  in  the  infeasible  region,  i.e.,  when  some 
of  the  constraints  are  not  satisfied. 

The  cam  base  diameter: 

In  the  cam-follower  mechanism,  it  is  generally 
always  desirable  to  have  the  smallest  cam[2].  Small  cams 
allow  us  to  build  a  compact  cam-follower  mechanism,  thus 
minimizing  problems  of  space  allocation.  Also,  the 
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follower  path  per  cycle  will  be  short  if  the  cam  size  is 
reduced,  which  in  turn  will  reduce  the  wear.  Vibration 
of  large  parts  is  always  undesirable,  and  the  use  of 
small  cams  also  reduces  the  chances  of  unbalancing  at 
high  speed.  Smaller  cams  are  also  cheaper  to  produce.  On 
the  other  hand,  if  the  cam  becomes  too  small  we  may 
encounter  large  pressure  angles,  high  contact  forces  and 
high  stresses[2]  .  We  may  also  run  into  problems  of 
undercutting . 

From  these  considerations,  it  is  clear  that  the  base 
diameter  of  the  cam  should  be  treated  as  a  design 
variable  in  the  optimization  problem.  The  cam  size,  as  a 
design  variable,  plays  an  important  role  in  both 
satisfaction  of  constraints  as  well  as  in  reduction  of 
cost. 

Orientation  of  lihe  cam  on  the  NC  table: 

In  the  previous  chapter  the  variable  step  generation 
technique  that  has  been  incorporated  into  this  work  was 
explained  .  Since  the  steps  generated  for  the  cutting 
tool  are  either  along  the  x-axis  or  the  y-axis,  these 
steps  are  very  strongly  dependent  on  the  orientation  of 
these  axes.   If  the  coordinate  axis  is  rotated  on  the  cam 
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surface,  a  different  set  of  steps  is  obtained  and  this 
has  a  direct  effect  on  the  manufacturing  time  .  For  this 
reason,  the  angle  representing  the  orientation  of  the 
coordinate  axes  is  treated  as  a  design  variable  .  This 
angle  of  orientation  of  the  axes  is  the  same  as  that  of 
the  orientation  of  the  blank  on  the  table  of  NC  machine. 
If  this  angle  is  initially  assumed  to  be  zero,  then  the 
orientation  of  the  coordinate  axis  will  be  such  that  the 
rise  polynomial  begins  from  the  positive  x-axis  .  Since 
the  step  size  and  number  of  turns  are  key  factors  in  the 
computation  of  the  cost  function,  this  design  variable 
plays  a  critical  role  in  reducing  the  manufacturing  cost. 
However,  it  does  not  affect  the  feasibility  of  the 
design . 

2)  COST 

The  minimization  of  a  cost  function  is  the  primary 
objective  of  optimization  .  As  mentioned  earlier,  the 
task  of  this  research  is  to  reduce  the  manufacturing 
time.  This  cost  function  utilized  here  can  be  broken  up 
into  two  parts  .  The  first  part  depends  on  the  tool 
diameter,  feed  rate,  total  distance  traversed  and  spindle 
speed,  whereas  the  second  part  is  comprised  of 
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the  time  lapse  per  turn  multiplied  by  the  total  number  of 
turns  in  the  tool  path  that  will  be  required  to  cut  the 
cam  .  The  details  of  the  estimation  procedure  for  the 
manufacturing  time  are  covered  in  Chapter  3 . 

3)  CONSTRAINTS 

Defining  a  proper  set  of  constraints  is  very 
critical  in  obtaining  a  usable  design  .  A  complete 
knowledge  and  understanding  of  the  problem  is  required 
before  the  constraints  are  identified.  Efforts  should  be 
made  to  review  all  performance  criteria.  Constraints  are 
the  only  restrictions  on  the  search  procedure  that 
prevent  the  search  from  converging  at  an  undesirable  or 
dangerous  solution. 

In  this  research,  constraints  are  imposed  to  prevent 
high  values  of  pressure  angle  and  contact  stress. 
Also,  undercutting  and  follower  separation  are  avoided  by 
imposing  suitable  constraints  .  Constraints  are  also 
imposed  to  keep  the  cam  profile  free  of  undesirable 
oscillations.  As  an  option, the  user  is  allowed  to  impose 
new  constraints  as  desired,  such  as  restricting  the 
maximum  values  of  jerk,  acceleration  and  velocity. 
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The   explicit   defining   inequalities   for   these 
constraints  are  the  following: 

1 )  Undercutting 

r(  0  )  -  1.5  *  r^   <=   0.  (46) 

where 

r(  0  )  is  the  instantaneous  radial  distance  to 

the  cam  profile 
8     is  the  angle  of  the  cam  rotation 
r-    is  the  radius  of  the  roller  follower 

2)  Elimination  of  oscillation  in  follower  motion 
For  rise 

-1  *  y'  <=  0.  (47) 

For  return 

y'  <=  0.  (48) 

where 

y  '   is  the  reduced  follower  velocity 


3 )  Pressure  angle 


♦  ^   -    $  ^   <=  0.  (49) 


where 

$  .   is  the  instantaneous  pressure  angle 

$  o   is  the  specified  maximum  pressure  angle 
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4)  Contact  Stress 


where 


o  ^/        a^^   -   1   <=  0.  (50) 


a .   is  the  instantaneous  contact  stress  on 


the  cam  contour 


a o   is  the  specified  design  stress 


5)  Follower  separation 


-1  *  F^  <=  0.  (51) 


where 

Fj   is  the  driving  force 

The  preceding  constraints  are  required  to  be 
satisfied  at  all  times  during  a  cycle  .  Checking  for 
constraint  violations  at  very  short  intervals  usually 
increases  the  computing  time  significantly;  at  the  same 
time  if  these  intervals  are  too  long,  then  it  is  always 
possible  that  the  optimization  process  might  miss  some 
constraint  violations. 

In  this  work  each  rise  and  return  cycle  is  divided 
into  15  equally  spaced  nodes  on  which  the  constraints  are 
evaluated  .  This  number,  however,  can  be  changed  for  a 
more  precise  check  on  constraint  violation.   The 
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constraint  functions  are  constant  all  through  the  dwell 
zones  .  For  this  reason  no  nodes  are  specified  in  the 
dwell  zone;  however,  the  constraints  are  evaluated  at 
each  terminal  of  each  zone  in  order  to  account  for  the 
behavior  of  constraints  in  the  dwell  zones. 

THE  OPTIMIZATION  STRATEGY 

The  optimization  strategy  is  based  on  the  concept  of 
an  exterior  penalty  function  method.  A  pseudo  -  objective 
function  "F  "  is  defined  as  follows: 

F   =  F  +  R  *  g  (52) 

where 

F    is  the  true  cost  function,  in  this  case  the 

cutting  time 
R   is  a  penalty  multiplier,  usually  specified 

to  be  of  the  order  of  10   or  10 
g   is  the  highest  constraint  violation 

The  minimum  cost  for  the  constrained  problem  is  then 
found  approximately  by  finding  the  unconstrained  minimum 
of  the  pseudo-objective  function[12]  .  Looking  at  the 
values  of  "g",  it  can  be  said  that  no  penalty  is  imposed 
if  all  constraints  are  satisfied;  however  if  one  or  more 
constraints  are  violated,  than  the  highest  violation  is 
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multiplied  by  "R"  and  added  as  a  penalty  to  the  cost 

* 
function  •   The  unconstrained  minimum  of  F   is  therefore 

likely  to  be  found  inside  the  feasible  region,  since  any 

constraint  violation  causes  a  sharp  increase  in  the 

pseudo-objective  function. 

The  cost  function  "F" ,  which  corresponds  to  the 
manufacturing  time,  is  discontinuous  with  respect  to  the 
design  variables .  This  is  true  because  the  manufacturing 
time  depends  on  the  number  of  turns  in  the  tool  path. 
Whenever  a  small  design  change  causes  a  change  in  the 
number  of  turns  of  the  tool  path,  the  cost  function  will 
undergo  a  step  change  .  Because  of  this  fact,  gradient 
based  methods  are  not  applicable  in  this  case.  The  method 
that  was  finally  used  to  perform  the  optimization  was  a 
directed  grid  search  technique[ 13 ]  which  does  not  require 
gradients  .  Although  it  is  a  relatively  crude  method,  it 
performed  adequately  well  in  the  current  work. 

OUTPUT 

Once  the  design  is  complete  and  manufacturing 
cost/time  is  optimized,  the  result  can  be  displayed  on 
the  screen  in  the  form  of  graphs  and  cam  profile  plots. 
The  final  task  of  the  process  is  the  generation  of  NC 
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codes,  which  can  then  be  transferred  to  metallic  or  paper 
tape.  The  tape  is  then  fed  to  the  NC  machine.  The  output 
also  provides  the  size  of  the  rectangular  blank  and 
orientation  of  the  cam  blank  on  the  table  of  the  NC 
machine . 
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Chapter  V 
NDMKRICAL  EXAMPLES 

The  techniques  developed  in  the  preceding  chapters 
was  implemented  in  a  computer  program.  In  order  to  verify 
the  approach  and  its  implementation,  several  numerical 
examples  were  solved,  some  of  which  are  presented  in  this 
chapter. 

For  these  examples,  machining  data  was  obtained  from 
Reference  [9]  and  the  physical  properties  were  taken  from 
Reference  [15].   The  machining  data  and  the  physical 
properties  were  kept  unchanged  in  all  the  examples. 
A  complete  set  of  these  constant  parameters  is  given 
below. 

Physical  Properties : 

Material: Hardened  Steel  SAE  1010 ( for ^both  cam  &  follower) 
Density  of  the  material  =0.28  lb/in 
Youngs  modulus  of  elasticity  =  29  x  10 
Poison's  Ratio  =  0.29 

Machining  Data: 

Tool  Material:  High  Speed  Steel 

Tool  diameter  =0.75  inches 

Number  of  teeth  =  4 

Cutting  speed  =  109.95  fpm 

Feed  per  tooth  =  0.0025  inch/tooth 

Spindle  speed  =  560  rpm_ 

Feed  rate  =  9.333  x  10~   inch/sec 
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Example  1 ; 

INPUT 

Design  Specifications ; 

Preload  =  320.0  lbs 

Maximum  Allowable  Stress  =  250000  psi 

Maximum  Pressure  Angle  =  12.0  degrees 

RPM  =  1000 

Mass  of  follower  mechanism  =  .02  Ibsm 

Lift  =  1 .20 

Cam  Thickness  =  0.25  inches 

Follower  diameter  =  0.85 

Follower  thickness  =  0.2  5 

Required  accuracy  =  l.OOOE-02 

Cycle  Characteristics ; 

Rise        Dwell        Return     Dwell 
0.0-180.0    180.0-185-0   185.0-350.0  350.0-360.0 
Boundary  Conditions 


'=0 


RISE 

At  z=0. 

y=o. 

y'=o, 

y"=o. 

At  z=l  . 

y=i. 

y'=o, 

y"'=o, 

RETURN 

At  2=0. 

y=i. 

y'=o. 

y'"=o. 

At  z=l 

y=o. 

y'=o. 

y"=o. 

Design 

Variables : 

'=0 


Number  of  design  variables  =  4 

Physical  significance  of  design  variables: 

B(l)  is  the  leading  coefficient  of  the  rise  polynomial. 

B(2)  is  the  leading  coefficient  of  the  return  polynomial 

B(3)  is  the  base  radius  of  the  cam. 

B(4)  is  the  orientation  of  the  cam. 
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Design  Variable  Initial  Value  Upper  Bound  Lower  Bound 

B(l)  4.0  180.         -180. 

B(2)  -0.072  180.  -180. 

B(3)  1.5  10.  1. 

B(4)  12.7  180.  -180. 


COMMENTS 

The  design  procedure  started  from  the  infeasible 
region  with  maximum  constraint  violation  of  5.99  55;  the 
pressure  angle  constraint  was  found  to  be  the  most 
violated.   When  the  design  stepped  into  the  feasible 
region  for  the  first  time  the  cost  function  was  found  to 
be  9  27.2  seconds. 


FINAL  RESULTS 


Design  Variables 
B(l)  =  3.996 
B(2)  =  2.0160 
B(3)  =  2.7398 
B(4)  =  18.684 


Coefficients  of  polynomials 

RISE  RETURN 

1  0.00  1.00 

2  0.00  0.00 

3  0.00  -4.495 

4  0.00  0.00 

5  5.005  9.959 

6  -1.400E-2  -5.919 

7  -7.987  -2.560 

8  3.996  2.016 
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Maximum  Values 

Force 
Velocity 
Pressure  Angle 
Contact  Stress  Sqr 
Positive  Accel 
Negative  Accel 
Positive  Jerk 
Negative  Jerk 

Highest  Constraint  = 

Cutting  Time  =  795.01  sec 

Number  of  turns  =  1298 

Tool  path  distance  =  25.7430  inch 

%  reduction  in  cost  from 

the  first  feasible  design  =  14.2587 


Magnitude 

Location  on  cam  profile 

(degrees ) 

.486E+03 

291.071 

.833E+02 

128.570 

.118E+02 

115.714 

.316E+11 

279.286 

.813E+04 

291.071 

■.160E+05 

180.000 

.202E+07 

232.143 

•.228E+07 

-10 

141.429 

■■   2.223368  * 

10  ^^ 

7t 


PesttTon      y 

V«loolty       y' 

Aoe«l«r        y" 

J-rK                 y'" 

L 

^jp^ 

^_y>- 

y,^^" 

(CW""^ 

\ 

\i 

RISE      CVC 

JLE 

^ 

J 

1 

Figure  12(a):    Example  1:  Results 
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Position    y 
V«loolty     y' 
Aoooi«r     y 
Jork            y'" 

/?K 

RETURN 

cvcle: 

L 

Figure  12(b):    Example  1:  Results 
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Example  2 ; 

INPUT 

Design  Specifications ; 

Preload  =  330.0  lbs 

Maximum  Allowable  Stress  =  250000  psi 

Maximum  Pressure  Angle  =  20.0  degrees 

RPM  =  1000 

Mass  of  follower  mechanism  =  .02  Ibsm 

Lift  =1.00  inch 

Cam  Thickness  =  0.25  inches 

Follower  diameter  =  0.75 

Follower  thickness  =  0.25 

Required  accuracy  =  l.OOOE-02 

Cycle  Characteristics ; 

Rise         Dwell       Return       Dwell 
0.0-150.0    150.0-160.0   160.0-330.0   330.0-360.0 
Boundary  Conditions 


RISE 

At  z=0. 

At  z=l  . 

y=o. 

y'=o. 

y"=o. 

y'"=0 

RETURN 

y=i. 

y'=o. 

y"=o. 

y"'=0 

At  z=0. 

At  z=l 

y=i. 

y'=o. 

y"=o, 

y   =0 

y=o. 

y'=o, 

y"=o. 

y   =0 

Design 

Variables : 

Number  of  design  variables  =  4 

Physical  significance  of  design  variables: 

B(l)  is  the  leading  coefficient  of  the  rise  polynomial. 

B(2)  is  the  leading  coefficient  of  the  return  polynomial 

B(3)  is  the  base  radius  of  the  cam. 

B(4)  is  the  orientation  of  the  cam. 
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Design  Variable   Initial  Value   Upper  Bound   Lower  Bound 


B(l) 
B(2) 
B(3) 
B(4) 


0.0 
0.0 
1.5 
0.0 


180. 

-180 

180. 

-180 

10. 

1 

180. 

-180 

COMMENTS 

The  design  procedure  started  within  the  feasible 
region  with  the  initial  cost  of  533.1  seconds. 


FINAL  RESULTS 


Design  Variables 

B(l)  =  0.000 

B(2)  =  0.000 

B(3)  =  1.45689 

B(4)  =  12.168 

Coefficients  of  polynomials 


RISE 

RETURN 

1 

0.00 

1.00 

2 

0.00 

0.00 

3 

0.00 

0.00 

4 

-2.168E-19 

-4.3368E-19 

5 

35.000 

-34.999 

6 

- 

-84.000 

83.999 

7 

70.000 

-69.999 

8 

- 

-20.000 

19.999 

9 

0.00 

0.00 

Maximum  Values 

Magnitude 

Location  on  cam  profi 

le 

(degrees ) 

Force 

.585E+03 

42.857 

Velocity 

.875E+02 

75.00 

Pressure 

Angle 

.198E+02 

64.286 

Contact  ! 

Stress 

Sqr.  .384E+11 

42.857 

Positive 

Accel 

.120E+05 

42.857 

Negative 

Accel 

-.120E+05 

107.143 

Positive 

Jerk 

.255E+07 

21.429 

Negative 

Jerk 

-.336E+07 

75.000 
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Highest  Constraint  =  3.55739019  x  10~°^ 

Cutting  Time  =  488.8689  sees 

Number  of  turns  =  800 

Tool  path  distance  =  15.7610  inch 

%  reduction  in  cost   =  8.45155 
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Peattton      y 
Valoolty       y 
Aoa«l«r        y*? 
U«rK                y'" 

^.A 

RISE      CYC 

JL-E 

i..e. 

Figure  13(a):    Example  2:  Results 
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Posftten      y 
V«loolty       y' 
Aao«l«r        y" 
J-rK                 y"* 

A 

RETURN 

CVCLE 



Figure  13(b):    Example  2:  Results 
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Example  3 ; 

INPUT 

Design  Specifications ; 

Preload  =  320.0  lbs 

Maximum  Allowable  Stress  =  250000  psi 

Maximum  Pressure  Angle  =  18.0  degrees 

RPM  =  1000 

Mass  of  follower  mechanism  =  .02  Ibsm 

Lift  =  1 .50 

Cam  Thickness  =  0.25  inches 

Follower  diameter  =  1.00 

Follower  thickness  =  0.2  5 

Required  accuracy  =  l.OOOE-02 

Cycle  Characteristics : 

Rise         Dwell        Return       Dwell 

0.0-190.0    190.0-190.0   190.0-360.0   360.0-360.0 

Boundary  Conditions : 


'=0 


RISE 

At  z=0. 

y=o. 

y'=o. 

y"=o. 

At  z=l . 

y=i, 

y'=o. 

y'"=o, 

RETURN 

At  z=0. 

y=i, 

y'=o. 

y'"=o. 

At  z=l 

y=o. 

y'=o. 

y"=o. 

Design 

Variables: 

'=0 


Number  of  design  variables  =  4 

Physical  significance  of  design  variables: 

B(l)  is  the  leading  coefficient  of  the  rise  polynomial. 

B(2)  is  the  leading  coefficient  of  the  return  polynomial 

B(3)  is  the  base  radius  of  the  cam. 

B(4)  is  the  orientation  of  the  cam. 
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Design  Variable   Initial  Value    Upper  Bound   Lower  Bound 


B(l) 
B(2) 
B(3) 
B(4) 


0.0 
0.0 
1.5 
0.0 


180. 

-180. 

180. 

-180. 

10. 

1. 

180. 

-180. 

COMMENTS 

The  design  procedure  started  from  the  infeasible 
region  with  maximum  constraint  violation  of  23  .75;  the 
force  constraint  was  found  to  be  the  most  violated.  When 
the  design  stepped  into  the  feasible  region  for  the  first 
time  the  cost  function  was  found  to  be  1388  seconds. 


FINAL  RESULTS 


Design  Variables 
B(l)  =  1.368 
B(2)  =  3.204 
B(3)  =  3.6548 
B(4)  =  28.872 


icients 

of  polynomials 

RISE 

RETURN 

1 

0.00 

1.00 

2 

0.00 

0.00 

3 

0.00 

-4.199 

4 

0.00 

0.00 

5 

8.29 

6.999 

6 

-9.212 

1.99E-: 

7 

0.554 

-7.015 

8 

7.368 

3.204 
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Maximum  Values 

Magnitude 

Location  on  cam  profile 
(degrees ) 

Force 

.651E+03 

299.286 

Velocity 

.940E+02 

135.715 

Pressure  Angle 

.104E+02 

122.143 

Contact  Stress  Sqr 

.292E+11 

299.286 

Positive  Accel 

.973E+04 

299.286 

Negative  Accel 

-.160E+05 

192.0 

Positive  Jerk 

.217E+07 

238.571 

Negative  Jerk 

-.204E+07 

149.256 

Highest  Constraint  =  2.63293682  10 
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Cutting  Time  =  1053.88  seconds 

Number  of  turns  =  1718 

Tool  path  distance  =  34.1495  inches 

%  reduction  in  cost  from 

the  first  feasible  design  =  24.129 
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1 

Poaftlon     y 
V«leelty      y' 
Aoe«l«r       y" 
J«rl<               y" 

RISE      CVC 

:i_E 

1.^; 

Figure  14(a):    Example  3:  Results 
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Pesttton   y 
V«leolty    y' 
Aooviar     y" 

/                          \v" 

RETURN 
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Figure  14(b):    Example  3:  Results 
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Chapter  6 


CONCLUSION 


The  motivation  for  the  work  presented  in  this  thesis 
was  the  perceived  need  for  new  techniques  in  the  area  of 
design  for  producibility  .  The  concept  of  accounting  for 
manufacturing  considerations  at  the  design  level  is  a 
very  significant  aspect  of  this  work  and  addresses  a 
major  weakness  of  Computer  Integrated  Manufacture  as 
implemented  in  industry  today. 

The  goal  set  for  this  work  was  to  achieve 
integration  of  design  and  manufacture  into  a  unified 
framework  and  to  account  for  manufacturing  considerations 
in  a  convenient,  general  manner  ,  A  mathematical  model, 
based  on  nonlinear  optimization  concepts  was  developed; 
this  model  makes  it  possible  to  handle  both  design  and 
manufacturing  considerations  simultaneously  .  Based  on 
this  model,  a  solution  technique  is  also  developed. 

The  use  of  nonlinear  optimization  techniques  offers 
distinct  advantages.  A  wide  variety  of  problems  can  be 
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easily  acconunodated  .  Both  design  and  manufacturing 
parameters  can  be  conveniently  included  in  the  problem 
statement  and  solution. 

The  solution  technique  for  the  nonlinear  optimiza- 
tion problem  is  a  systematic  iterative  process.  Starting 
from  any  initial  design,  the  technique  produces  a 
sequence  of  improved  designs  until  the  optimal  solution 
is  reached  .  Furthermore,  the  technique  can  be  easily 
implemented  in  a  computer  code. 

Another  major  advantage  of  this  approach  is  that  the 
design  process  becomes  highly  automated.  Once  the  initial 
values  of  the  design  variables  are  specified  and  the  cost 
and  constraint  functions  are  defined,  the  process  finds 
the  optimal  design  with  no  further  user  interaction.  This 
not  only  cuts  down  the  amount  of  designer  time  that  is 
required  but  also  reduces  the  level  of  skill  required  to 
produce  the  design. 

The  model  also  guarantees  an  improvement  in  the 
initial  design  .  With  any  initial  design  the  solution 
first  proceeds  in  the  direction  of  constraint  correction 
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and  then  in  the  direction  of  cost  reduction.  It  provides 
at  least  as  good  a  design  as  that  which  would  be  produced 
by  a  human  designer  using  conventional  design  techniques. 

The  implementation  of  the  proposed  solution 
technique  proved  to  be  very  reliable  and  computationaly 
viable  as  can  be  seen  from  the  numerical  examples  .  The 
examples  also  point  out  the  generality  of  formulation. 
The  model  can  handle  any  number  of  boundary  conditions  as 
well  as  any  specified  values  for  the  lift,  rise/return 
zone  widths,  maximum  pressure  angle,  maximum  stress,  etc. 

Although  the  results  were  very  encouraging,  there 
are  some  serious  limitations  to  this  work: 

1)  Manufacturing  considerations  have  not  been 
treated  in  great  detail . 

2)  Only  plate  cams  with  roller  followers  are 
considered . 

3)  Only  single  rise  and  single  return  cam  profiles 
are  considered. 

4)  The  numerical  technique  that  has  been  used  is 
somewhat  crude,  even  though  it  has  been 
effective. 
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Future  work  is  required  to  address  the  above 
limitations  .  The  model  should  be  expanded  so  that 
manufacturing  parameters  such  as  cutting  speed,  tool 
diameter,  feed,  etc  .,  can  also  be  treated  as  design 
variables.  Tolerances  on  the  machined  surface  should  be 
given  more  consideration,  and  studies  should  be  done  on 
how  these  tolerances  relate  to  the  manufacturing 
parameters.  The  use  of  piecewise  polynomial  definitions 
(e  .g  .,  B-Splines)  instead  of  the  single  polynomial 
approximation  used  in  this  thesis,  will  make  it 
possible  to  design  cams  that  can  satisfy  more  rigorous 
performance  constraints  and  provide  greater  reduction 
in  the  cost. 

The  model  can  also  be  extended  to  include  other 
type  of  cams/followers  and  to  account  for  cases  with 
multiple  rises  and  multiple  returns.  Finally,  there  is 
room  for  improvement  in  the  numerical  methods  that  have 
been  used  .  These  extensions  will  make  the  model  more 
comprehensive  and  greatly  expand  the  range  of 
applications  for  which  it  can  be  used. 
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ABSTRACT 


The  objective  of  this  work  is  to  integrate  the 
design  and  manufacture  of  Plate  Cams  through  the  use  of 
non-linear  optimization  techniques  .  The  underlying 
concept  is  design  for  producibility .  A  constrained  non- 
linear minimization  problem  is  set  up  in  which  the  cost 
function  is  taken  to  be  the  manufacturing  time  required 
to  produce  the  cam .  Constraints  are  imposed  to  ensure 
that  the  design  is  acceptable  with  respect  to  criteria 
such  as  pressure  angle,  undercutting,  stresses 
developed,  etc  .  The  design  variables  are  the  base 
radius  of  the  cam,  the  orientation  of  the  cam  blank  on 
the  NC  machine  table,  and  the  leading  coefficients  of 
the  polynomials  which  describe  the  rise  and  return 
cycle  of  the  cam.  A  grid  search  method  is  used  to  find 
the  constraint  optimum  iteratively  .  This  iterative 
scheme  was  successfully  implemented  in  a  FORTRAN 
computer  program  .  At  each  iteration  the  tool  path 
corresponding  to  the  current  cam  design  is 
automatically  generated  .  The  cost  and  constraint 
function  are  also  evaluated  automatically  .  After  the 
optimum  is  found,  the  results  are  displayed  graphically 
and  in  printed  form,  and  NC  code  required  for  producing 
the  cam  is  also  generated  •  Examples  are  presented  to 
demonstrate  the  effectiveness  of  the  approach  in 
automatically  designing  plate  cams  that  are  optimal 
with  respect  to  design  and  manufacture. 


